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Abstract 

We introduce an integrable time-discretized version of the classical Calogero- 
Moser model, which goes to the original model in a continuum limit. This 
discrete model is obtained from pole solutions of a discretized version of the 
Kadomtsev-Petviashvili equation, leading to a finite-dimensional symplectic 
mapping. Lax pair, symplectic structure and sufficient set of invariants of the 
discrete Calogero-Moser model are constructed. The classical r-matrix is the 
same as for the continuum model. 
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1. Introduction. The original Calogero-Sutherland-Moser model, JTJ, 0, [|, is a 
one-dimensional many-body system with pairwise inverse square interactions, and 
it is an integrable model at both classical and quantum levels, cf. the two review 
papers, for a comprehensive account of both cases. (We will refer to the classical 
model as the Calogero-Moser (CM) model). This long-range interaction model has 
aroused renewed interest most recently in view of its close relations to the Haldane- 
Shastry SU(N) chain that has been studied in ||. Intriguingly, there is strong 
evidence that this integrable model plays an important role in understanding the 
universal behaviour in quantum chaos and mesoscopic physics, ||. Furthermore, 
rich algebraic structures, such as a dynamical r-matrix, have been found in this 
model in 0, cf. also [§]■ 

In this letter we introduce a time-discretized version of the classical CM model. 
Discrete-time systems are a subject of intensive research recently, (cf. e.g. || 10 



for reviews), and they are, in our belief, important in understanding the true nature 
of integrability. Most known integrable systems, such as integrable partial differen- 
tial equations (Korteweg-de-Vries equation, Kadomtsev-Petviashvili equation, etc.) 
or ordinary differential equations (the Painleve transcendents and elliptic function 
equations), possess, as it has turned out to be the case, one or more integrable 
discrete counterparts which lead to the original model in a well-defined continuum 
limit. Finding such integrable discrete versions, however, is in general a highly 
non-trivial undertaking. In the present paper, we propose a new discrete integrable 
system, namely an exact discrete-time version of the Calogero-Moser system. 

The method used here to obtain the discrete-time CM model is based on the 



observation of [JIT], cf. also [12|, that the dynamics of the poles of special solutions 
of integrable nonlinear evolution equations is connected with integrable systems of 
particles on the line. The connection between the pole solutions of the (continuum) 
Kadomtsev-Petviashvili (KP) equation and the CM system was found by Krichever 
in [13], cf. also Here we will perform a similar construction for the discrete 
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case, and we will show that the pole solutions of a semi-discretized version of the 
KP equation is connected with a time-discretized version of the CM model. 

2. Discrete-time CM model. The semi-discretized version of the KP equation 
that we will use reads 

(p — q + u — u) y = (p — q + u — u){u + u — u — u) , (1) 

where p and q are two (lattice) parameters, u is the (classical) field, the ~ denotes 
the discrete time-shift corresponding to a translation in the "time" direction while 
the " denotes a shift or translation in the "spatial" direction. In (JTJ) the time and 
spatial variables are discrete while the third variable y is continuous. Eq. ([[]) is 
obtained from the fully-discretized version of the KP equation of [15| by letting one 
of three lattice parameters tend to zero. 

It can be easily proved that (|T|) is the compatibility condition of the following 
two equations 

= (j) y + (p + u - u)(j), (2a) 
= 4>y + (q + u - u)(j), (2b) 

which is the Lax representation of ([!]), and the starting point for the investigation 
of the integrability characteristics of the semi-discrete KP. 

Direct calculation shows that 

V, — — with x = — I h y (3) 

x p q 

is a (pole) solution to the semi-discretized version of KP equation ([I]). Here n and m 
are two integers, and (n, m) are considered as the coordinates of a two-dimensional 
lattice and 

_ n + 1 m n m + 1 

x = 1 h y, x= — I h y. (4) 

p q P q 

If we substitute this special solution @ into the it's in (§), then we find that 

<p= (l-^-)(p + k) n (q + k) m exp(ky) (5) 

rCX 
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satifies (0), where k is a spectral parameter. 

Enlightened by the above simple exercise, we now suppose [] that 

N 1 

u = , (6a) 

1 N b 

-^—)(P + k) n (q + k) m exp(ky), (6b) 

where and 6j are independent of y (but they depend on the time variable), and 
we are to find the conditions these Xi and 6j should satisfy such that equation (|2a|) 
is valid. Substituting (|6p into (|a|) and equating to zero the coefficients of (y — Xi)' 1 



and (y — xi) , we obtain the following equations: 



N b N b +b- 

■ ™ /yi . rp . * " /y> . /y> . 



- N b- 

(p + k)bi = k-J2=^, (7b) 

j=l x i x j 

for alH = 1, 2, N, where fej and denotes the discrete time-shift as stated before. 
If these conditions (0) are satisfied, then the u and given by (H) satisfy (^a|). 

Now, by introducing the vectors B = (pi, b 2 , &tv) t and i? = (1, 1, 1, 1) T and 
the matrices 

/ 



N I N 



E = " E % C 1 " , (8a) 

C * rp . rp ' * rp . v< , I J rp . rp , •* ' v ' 



I . 1 



MV=— -^.--—-(l-dy), (8b) 

we can rewrite (^|) in the following form 

(p + = kE + LB, (9a) 

(p + fc)5 = kE + MB, (9b) 



1 We note here that if we take the more restrictive Ansatz cj> = jgQ^Li y - x - ) ex P(ky) instead of 
that given in (Eh we still get the same result, i.e., eqs. ( |TT| ) and ( |T^ ) below. 
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and the compatibility of (|S|) leads to the equation 

(LM - ML)B + k(L - M)E = 0. (10) 

([T0"| ) is a discrete non-homogeneous Lax's equation. It can be readily checked that 
the resulting equations of (|10|), i.e. 

LM = ML (11) 

and 

(L-M)E = 0, (12) 

are consistent and give the same discrete equations of motion of a N— particle sys- 
tem: 




i = l,2,...,JV, (13) 

where Xi denotes the discrete time-shift in the opposite direction to the one of x^. 
We will call the model, for which eq. ([IB]) are the equations of motion, the discrete- 
time CM model. We will show below that in a continuum limit these equations go 
to that of the original CM model. 

The Lax pair for this discrete-time CM model are given by L and M of From 
eq. (|10"D , it can be readily seen that 

h = I k (14) 

for any k — 1, 2, where 

I k = Tr(L k ), (15) 

leading to a sufficient number of invariants (or conservation laws) of the discrete- 
time flow given by (|T3|). 
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3. Symplectic Structure. In order to establish the exact integrability of the dis- 



crete CM model fll3"D , we need first to establish an appropriate symplectic structure 
for the N-particle system. The discrete-time flow will then have the interpretation of 
the iterate of a canonical transformation with respect to that symplectic structure. 
To get the generating function of this canonical transformation, we will follow the 



point of view of ref . [[LJJ . We start by noting that eq. (|T3"D can actually be obtained 
from the variation of a discrete action, given by 

\ 

(16) 



s = j:a 



X,X) 



N N 

E | E hg\xi-xj\ - lo£ 



I X 2 X j I 



J 



in which the sum over n denotes the sum over all discrete-time iterates. The discrete 
Euler-Lagrange equations 



d£ + dC_ 

dx; dxl 



0. 



l,2,..,iV, 



(17) 



yield the equations of motion ([13|) . 



The form of the 'kinetic' term in our action (|1^) invites to perform a Legendre 
transformation in the following form 



N 

E 



H(p, X) = E Pij( X i ~ X j) ~ £( X > X ) 



from which we obtain 



and 



A? 



E^y 



dC 



i=l 



m __ 



dpi 



(18) 
(19) 

(20a) 
(20b) 



Eqs. fl20a|) and ( p0b| ) can be interpreted as discrete Hamilton equations, but, of 
course, TL is not a Hamiltonian in the usual sense of the word. H is the generating 
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function of the canonical transformation x i— > x, p h- > p. This transformation, as a 
consequence of eqs. (p0|) , will leave the following symplectic form invariant: 

N 

Q = dpij A dxj, (21) 

for which we have Q = Q. 

Eqs. (|i8|) -(pTl) are still general. Coming back to the special case of the action 



([16]) we will find from (0) that 



p- = — L_ (22) 



and from (|1|) we find 7i: 



N N 

tt(P,x) = l °g\Pij\ + log | ^ | . (23) 



i,j = l »>J'=1 



The symplectic structure ( plf ) leads to the following Poisson brackets 



N 

YliPii^k) = $jk- (24) 
1=1 

4. Classical r-matrix. In order to prove the complete integrability of the discrete- 
time CM model, we now construct its r-matrix structure. From the above discussion 
it follows that we can make the following choice of canonical variables (pi, Xi) of the 
ith. particle 

N N \ 

Pi = HPii + Y, > ( 25 ) 

7 -=i -^i 

leading to the standard Poisson brackets 



j=l j=l i 



{pi,x j } = 5 ij , {pi,pj} = {xi,Xj} = 0. (26) 



In terms of these canonical variables, the Lax matrix L can be written as: 

v e 



Laptev - 52 -22-, (27) 



i=l i,i=l X i X 3 



7 



where the elements of the matrix are defined as (e^-)jy = 8^8^. Eq. fl2"7|) is the 
usual L-matrix for the CM model, and thus we can immediately use the result of 
in order to establish the involutivity of the invariants flI5|). In fact, using (26) and 
the expression for L, we can calculate the fundamental Poisson bracket structure in 
terms of the matrices L, leading to the well-known result 

{LfL} — [n 2 , L®1\ — [r 2 i, 1 ® L], (28) 

where ® is the tensor product, [•, •] denotes the usual commutator and 

N N 

{LfL}=Y^ E i L ih L ki}eij ® e H . 
i ,j=i fc,/=i 

The r-matrix in (^S|) is given by, [|7|], 

w i 1*1 

r i2 = E ^ _ ^ % ® e i< + 9 E ~ _ z e a ® ( e y - e ^)> ( 29 ) 

obeying r 21 = Pr 12 P, where P is the permutation matrix: Px £g> yP — y®x. 

Thus, in going from the continuous to the discrete CM model, both the L- 
operator as well as the r-matrix remain the same. What changes is the M-matrix, 
which will now depend on the original momentum variables pij of eq. (p2|). 

The involutivity of the invariants 

{I k J l } = {Tr(L k ),Tr(L l )} = for all k, I = 1, 2, .... , (30) 



that follows as a consequence of the r-matrix structure (28), will lead to the in- 
tegrability of the discrete-time model by an argument presented elegantly in 
forming the discrete counterpart of the Arnol'd-Liouville theorem. Therefore, we 
have proved the complete integrability of the discrete-time CM model. 



5. Continuum Limit. We now show that in a continuum limit (|13|) goes to 
that of the original CM model. In fact, the considerations above show that the 
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invariants and the canonical variables remain the same on the continuous- as well 
as the discrete- time level. Thus, taking the invariants as Hamiltonians, we get 
a hierarchy of continuous flows interpolating the discrete-time flow, the one of order 
k = 2 corresponding to the original model. In order to perform the continuum limit, 
we first set 

Xi = Zi + na, i — l,2,...,N, (31) 
where a is a small (constant) parameter and 

xl=z-i + {n + l)a, i = l,2,...,N. (32) 

Then, in the continuum limit, we write 

zi = Zi + eii + \/2(?Zi + (33) 

for all i = 1,2,..,N, where ij = 4zj and e is the time-step parameter which, we 
suppose, is in the order of 0(a 2 ). Substituting fl3~i"])- (|33|) into (|13|), we get as the 
leading order term of fll3h, i.e. 

N 1 

Zl = -2gJ2 ( _ , 3 , i = l,2,...,N, (34) 

j=i [Zi Zj) 



where g = a A /e 2 . It is clear that (34) are exactly the equations of motion of the 
(continuous) CM model. It is interesting to note that the coupling constant of the 
continuous model arises from the discrete-time step. 



6. Disscusion. In this letter, we have constructed a discrete-time CM model. 
This discrete model is also integrable and is the iterate of a canonical transforma- 
tion generated by the discrete 'Hamiltonian' (^). The original CM model is a limit 
of this discrete-time CM model, and in this limit the coupling constant for the long- 
range interaction term is encoded in the discrete-time step parameter. The present 
result invites a number of interesting problems to be studied. First of all, one should 
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address the generalization of these results to the more generic elliptic potentials, cf. 



13, M. Secondly, one may ask whether one can find time-discretizations of the rela- 



tivistic version of the CM model, [T7||, cf. also (T^|. In fact, Suris, in |19[], has found 
an interesting connection between the discrete-time Toda model and its relativis- 
tic version. Furthermore, an intriguing connection exists between the sine-Gordon 



soliton solutions and the relativistic CM model, cf . [ 2(| . This could lead to an alter- 
native way to discretize the model. Thirdly, for the discrete-time model, following 
the similarities with the structure of the integrable quantum mappings studied in 
PH , one should investigate not only the L-part of the Lax pair, but also take the 
M-part under consideration in the r-matrix structure. This has been done for the 



integration of mappings of KdV type in p2| . Finally, these investigations should 
also be pursued on the quantum level. Although, the discrete-time model has an 
obvious quantum counterpart, in much the same way as the continuum model, the 



work of [21[ on quantum mappings, as well as the work |23| on the quantum CM 
model, indicate that some important modification (e.g. with respect to the con- 
struction of exact quantum invariants of the discrete-time flow) might be expected. 
All these problems are being investigated and will be discussed in detail in future 
publications. 
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